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Abstract: The Sornette-lde differential equation of herding and rational 
trader behaviour together with very small random noise is shown to lead to 
crashes or bubbles where the price change goes to infinity after an unpre- 
dictable time. About 100 time steps before this singularity, a few predictable 
roughly log-periodic oscillations are seen. 

Bubbles and crashes are a property of market prices since centuries. Some 
may arise from external perturbations like changes of the interest rate by the 
central bank, government decisions on war and peace, or for stocks of single 
companies the introduction of new (un) successful products. Some, like the 
October 1987 crash on Wall Street, may arise from intrinsic market mecha- 
nisms [1]. While we cannot make a clear distinction between external and 
intrinsic reasons for a crash in reality, we can do so at least in computer sim- 
ulations. There a single external event can be put in explicitly into the pro- 
gram, the multitude of individual decisions can be approximated by a small 
random noise while the general rational as well as psychological ( "herding" ) 
behaviour of investors and traders can be approximated by deterministic al- 
gorithms. We show here that in the Sornette-lde [2,3] approach this small 
noise can escalate to a singularity (bubble or crash) where the price change 
goes to ±cx). 

The Sornette-lde approach [2] used two different contributions. One is 
the nonlinear fundamentalist assumption of [3] that traders buy (sell) if the 
prices are low (high). Thus a positive difference x between the actual and the 
perceived fundamental price encourages selling and thus causes a downwards 
tendency in v = dx/dt such that dv/dt is proportional to —a;". This effect 
stabilizes prices. But a rising price, f > 0, creates the hope or illusion that the 
price will increase further and thus will encourage further buying, with dv / dt 
varying as v"^; this second term alone leads to a divergence v oc (tc— 1)^^^*^™~^^ 
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at some critical time t = tc (bubble or crash). The combination of both 
terms leads to roughly [4] log-periodic oscillations oc sin [const log(tf. — t)] 
preceding the singularity, which may be used to make profit or to prevent 
this singularity. 

One sample, m=3, n=5, x multiplied with 1 + 0.0001 * noise, -1 < noise < +1 
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Figure 1: Long-time behaviour of one example. Slow weak oscillations are 
followed by strong rapid oscillations. 

The deterministic nonlinear differential equation for the "return" x{t) (the 
logarithm of the current price P to the initial "fundamental" or "equilibrium" 
price) and its velocity (short-time return) v — dx/dt is 

dv/dt = CmV"'' — CnX^ 

if m and n are odd integers. The damped harmonic oscillator has m — n — 

1, Cm < 0, c„ > 0. Suitable absolute values are needed if m and n are 
not odd integers [2]. Noise can be introduced here in various ways: Chang 
et al [5] added a small random value to a;; then after a long time regular 
oscillations appeared, of increasing strength and decreasing period, until x 
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Three examples, m=3, n=5, return multiplied with 1 +0.0001 *noise, -1 < noise < +1 
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Figure 2: Behaviour shortly before the crash of the same sample as in 
Fig.l plus two other examples differing only by the random number seed. 
Only for a few oscillations the three curves are reasonable synchronized; for 
longer times before the crash some show a maximum where the others show 
a minimum. 

diverged. While this behaviour is what we want, the assumption about the 
noise has little justification except that it is close to Langcvin equations. 

We instead first made vary randomly between and +1, but then the 
curves were too regular, and different initial seeds for the random number 
generator gave nearly the same curves. This multiplicative assumption was 
supposed to simulate the volatile psychology of the traders who follow more 
or less the current trend. When we applied the random noise to c„ instead of 
Cm this noise infiuenced the prices more strongly but we got oscillations right 
from the beginning while we want them to emerge only later. Also applying 
noise to v did not give satisfactory results. 

To get better results, in our second method we made a multiplication for 
X itself, not for a contribution to its time derivatives. Thus the x resulting 
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Average over 5000 crashes 

4000 r- 
3000 - 
2000 - 
1000 - 



-1000 - 
-2000 - 
-3000 - 

-100 -80 -60 -40 -20 

time before crash 

Figure 3: Average over many crashes, for the same parameters as in Figs. 1,2. 

from the discretized differential equation was multiphed by a factor which 
differed from unity by a small amount between -0.0001 and +0.0001. Since 
we use multiplication instead of addition, we cannot start with x and v both 
zero, and thus took x = 0, v = 0.0001 initially, to create some fluctuations. 
We thus assume that traders decide in a deterministic way on fundamental 
value through Cnx"' and on herding through Cmv'^, but finally some small 
randomness also affects the market. 

Figure 1 shows a run over all nearly 5000 time steps. Fig. 2 is the same run 
over the last 200 time steps before the crash, together with two other runs 
using the same parameters but different random numbers. We see that first 
the prices fluctuate in a rather random way about the fundamental value. In 
other runs, the noise can make the deviations larger leading to a few nearly 
regular oscillations which die down again. Finally, again such oscillations 
occur, but now they become stronger, faster and finally lead to a crash at 
time = 0. Averaging over thousands of samples gives smooth oscillations 
before the crash. Fig. 3. If we increase the noise, the times needed for a crash 
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Histogram for lifetimes of market until crash; noise = 0.01 , 0.001 , 0.0001 from left to right 
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Figure 4: Distribution of the times after initialization needed for a crash 
to appear. A parabola in this log-log plot means a log-normal distribution. 
The bin sizes increase by a factor 2. Besides our normal noise 0.0001 on the 
right, also ten and hundred times stronger noise is shown (center and left 
curves.) 



to build up get shorter and there are less oscillations; these times are roughly 
log-normally distributed, Fig.4. 

Mathematically it is easy to understand why this second method works 
better than our first method. In the first method, the noise affects the time 
derivative only, meaning it changes the curvature or slope of x{t) only. Our 
second method affects x directly. 

Without any noise, the origin x — v — ^ oi the [y^ x) phase space plays 
a special role as the unstable fixed point around which spiral structures of 
trajectories are organized. It corresponds to the case of fundamental prices; 
there is no trend and the market does not know which direction to take. 

Our parameters for this simulation were m = 3, n = 5, = = 1. We 
used the leap frog method of molecular dynamics with 1000 iterations per 
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unit time. Thus first we calculate from the above equation the acceleration 
a = (Px/dt^, from a the change 0.001a in the velocity v = dx/dt, and finally 
the new x is obtained by first adding O.OOlf to it using the new velocity, and 
then multiplying x by 1 + e, with e selected randomly between -0.0001 and 
+ 0.0001. 

The "noise" e symbolizes all the new information coming in every day (or 
every simulation interval). We also tried to simulate insider trading: Half of 
the market movement uses not today's noise but the one of tomorrow. Not 
much is changed since the equations describe the overall market behaviour, 
not the profits of some at the expense of others. 

Of course, this method is similar to that of Chang et al; the main dif- 
ference is that the present noise in x is proportional to x while for Chang 
et al it is independent of x. An additional linear friction term, dv/dt = 
CmV^ — CnX"' — CfV, restricts the fiuctuations x mostly to one sign, except for 
very small friction coefficients c/. 
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